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Abstract 

We study the statistical mechanics of a multicomponent two-dimensional Coulomb 
gas which lives on a finite surface without boundaries. We formulate the Debye- 
Hiickel theory for such systems, which describes the low-coupling regime. There 
are several problems, which we address, to properly formulate the Debye-Hiickel 
theory. These problems are related to the fact that the electric potential of a single 
charge cannot be defined on a finite surface without boundaries. One can only define 
properly the Coulomb potential created by a globally neutral system of charges. 
As an application of our formulation, we study, in the Debye-Hiickel regime, the 
thermodynamics of a Coulomb gas living on a sphere of radius R. We find, in this 
example, that the grand potential (times the inverse temperature) has a universal 
finite-size correction (1/3) In R. We show that this result is more general: for any 
arbitrary finite geometry without boundaries, the grand potential has a finite-size 
correction (x/6) InR, with \ the Euler characteristic of the surface and R 2 its area. 
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1 Motivation 



In this paper, we study two-dimensional Coulomb systems, for instance plas- 
mas or electrolytes, which live on a finite surface without boundaries. The sim- 
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plest example of such geometry is the sphere. This kind of geometry has been 
used in numerical simulations of charged systems [1,2,3,4] as an alternative to 
the Ewald method, since it avoids the problem of the boundary conditions. 
There has been also several theoretical studies of Coulomb systems on the 
sphere, in particular, the statistical mechanics of both the one-component 
plasma and the charge symmetric two-component plasma, can be exactly 
solved for a special value of the Coulomb coupling [5,6,7]. The purpose of this 
work is to formulate the Debye-Hiickel theory for a generic multi-component 
Coulomb system living on a surface without boundaries. The Debye-Huckel 
theory describes the low-coupling regime of the system. 

There are several difficulties which make this problem non-trivial. First, it is 
not possible to define the Coulomb potential as the inverse of the Laplacian 
operator, since on a surface without boundaries, the Laplacian has no inverse. 
Physically, this means that the electric potential of a single charge, on the 
sphere, does not exist. One can only properly define the electric potential 
created by a globally neutral configuration of charges. Thus, the partition 
function of such system should be restricted to neutral configurations. This 
problem is not present on an infinite surface, or a surface with boundaries, 
where any excess charge can go to infinity or to the boundaries, and does not 
affect the bulk thermodynamics of the system. 

There are several equivalent ways to formulate the Debye-Huckel theory. Re- 
cently, the author and a collaborator proposed one, which is particularly ap- 
propriate to study confined Coulomb systems [8,9]. This method is based 
on the Hubbard-Stratonovich transformation, also known as the sine-Gordon 
transformation, when applied to Coulomb systems. However, for a Coulomb 
system living in a surface without boundary, this method is not directly ap- 
plicable, because: 1) to use the sine-Gordon transformation one should con- 
sider all possible configurations of the system, including the globally charged 
ones, 2) the Laplacian, which in a flat geometry appears as the inverse of the 
Coulomb potential, has no inverse on a surface without boundaries. 

The outline of this work is as follows. In the next section, we will show how to 
overcome the difficulties mentioned above and how the approach of Refs. [8,9] 
can be extended to geometries without boundaries. For simplicity we will 
consider the case of the sphere, however the method presented in Sec. 2 is 
general enough to be applicable to any surface without boundaries. In Sec. 3, 
using the results of Sec. 2, we will explicitly compute the grand potential and 
other thermodynamics functions of a Coulomb system living on a sphere in the 
Debye-Huckel regime. We will also compute the finite-size expansion of the 
grand potential, and show the existence of a finite-size correction (1/3) In R 
to the grand potential (times the reduced inverse temperature), where R is 
the radius of the sphere. Finally in Sec. 4, with the aid of some results from 
Ref. [9], we will show that for a general geometry without boundaries, the 
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grand potential has a logarithmic finite-size correction (x/Q)lnR where x is 
the Euler characteristic of the manifold where the system lives and R 2 is 
proportional to the area of the manifold. 



2 Modified sine-Gordon transformation 



2. 1 Coulomb potential and energy of a pair of pseudocharges 



The Coulomb potential is usually defined as the solution of Poisson equation 

Av(r,r') = -2tt6(t,t') (2.1) 

with appropriate boundary conditions. For a flat plane geometry, the potential 
is 

v°(r,r') = - In ^ ~ ^ (2.2) 

which satisfies (2.1) and the boundary condition Vi>°(r, r') — > as |r — r'| — > 
oo. In equation (2.2), L is an arbitrary constant which fixes the zero of the 
electric potential. 

On a sphere, of radius R, and more generally on any finite surface without 
boundaries, it is not possible to define the Coulomb potential from Poisson 
equation (2.1). This equation has no solution on the sphere: the electric po- 
tential of a single charge cannot be defined. Instead, following [3], we should 
consider that the system is composed of pseudocharges. A pseudocharge is the 
combination a unit point charge and a uniform background of opposite charge 
spread over the sphere. The electric potential v created by a unit pseudocharge 
satisfies 

Av(Q, SV) = -2vr (,5(fi, ft) - — ) . (2.3) 

with Cl = (9, ip) the spherical coordinates of source and Q' = (9', ip') the 
location where the potential is computed, and 8(0,, Q') = 5(cos 9 — cos 9')5(p — 
(f') I R 2 is the Dirac distribution on the sphere of radius R. 

Decomposing in spherical harmonics Y^ m , one can write down the solution of 
the previous equation as 

oo £ i 

v(Q;n')=2nY: £ ——-Y, m {9'^')Y, m (9, V ) + V (2.4) 
e=i m =-e l \ l + l ) 

with Vq an arbitrary constant. Now, using [10] 
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where Pe(x) are the Legendre polynomials of order £, and d is the angle be- 
tween Q and Q', we can write the Coulomb potential (2.4) as 

u(0,p;0>) = -lnsin^-i + Vo- (2.6) 

The potential energy Eij of two pseudocharges q^ and qj located at fli and Qj 
can be decomposed as = Vij + Si + £j, with the interaction energy and 
£j (£j) the self-energy of the pseudocharge (^). The interaction energy is 

^. = ft y V (Q., f2)pj(f2) R 2 dQl (2.7) 

where Pj(fi) = qjS(Q, Qj)—qj/ (AnR 2 ) is the charge density of the pseudocharge 
qj. We find 

^ = -^(msin^-^ (2.8) 

with the angle between Jlj and f^- measured from the center of the sphere. 
Notice that the interaction energy does not depend on the arbitrary constant 
V . 

Following [3], the self-energy £j of a pseudocharge can be computed by in- 
tegrating the square of the electric field E(f2) = — Vt>(f2,f2j) created by the 
pseudocharge on the whole sphere. However this self-energy is infinite for a 
point charge. Instead we can replace, initially, the point charge by a small disk 
of radius a with the charge % spread over its perimeter and we compute 

e t = ^ / \E(n)\ 2 d& - el . (2.9) 

47T J 

We have subtracted ef, the self-energy of the charge alone, without its 
neutralizing background. This self-energy ef is not properly defined on the 
sphere, since it correspond to a non-neutral charge configuration. As in Ref. [3] 
we adopt the prescription of taking ef as the self-energy of the charge q^ on a 
flat surface 

^ = -fm^ (2.10) 

In the limit of a point charge, a — > 0, the self-energy of the pseudocharge (2.9) 
is finite 

£ ,.|(_ 1 + 21n ^). ( , n) 

Let us now study a Coulomb system on the sphere composed of several species 
of pseudocharges q a , a — 1, • • • , s. Each species a has N a pseudocharges. The 
position of the i-th pseudocharge of the species a will be denoted by £l ai . The 
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total potential energy of the system is 

1 N <* N -> _ N 

H = o E E E 'w^fc, n™) + EE £ « (2.12) 

^ a, 7 j=l jr'=l a i=l 

with 

Uc^iAj) = -lnsin(^/2) - 1/2 (2.13) 

and dij the angle from the center of the sphere between Vt a i and fi 7 j. The 
prime in the summation in (2.12) means that the term when a = 7 and % — j 
must be omitted. 

We shall call "neutral" any configuration of pseudocharges satisfying J2 a Q a N a = 
0. For a neutral configuration, the backgrounds of the pseudocharges cancel 
each other, and the system is equivalent to a system of charged particles 
without neutralizing backgrounds. For those neutral configurations, we have 

(£ Q £ l= ?i Qa) 2 = = E a £& q 2 a + E a , 7 £*=" Efji W Th en, the potential 
energy can be put in the form 

^neutral = ~ \ E E '?a?7 ln ( ~T sin ^ ) " ( 2 ' 14 ) 

In the flat limit R — > oo, keeping = iMjj finite, the above expression 
reduces to iffl at = —(1/2) £ ai7 g a g 7 ln(ry/L), which is the hamiltonian of 
a Coulomb system in the flat geometry. Thus, the prescription (2.11) for the 
self-energy is a reasonable one, if one wishes to recover the results for the flat 
geometry in the limit R — > oo [4]. Notice, however, that the thermodynamics of 
the Coulomb system on the flat geometry can only be recovered if we consider 
only neutral configurations from the start on the sphere. 

2.2 Sine-Gordon transformation 

We shall work in the grand canonical ensemble. The fugacity of the species a 
will be denoted by ( a . We define as usual (3 = 1/(/cbT), were T is the absolute 
temperature and ks is the Boltzmann constant. The grand partition function 
reads 

oo oo fNi fN s r r r N a 

s = E • • • E fr 1 ^ / • • • / ^ H n n (2- 15 ) 

Af 1= N s =0 l\i....l\ s .J J Q =li=l 

The Kronecker symbol 5y^ N ensures that we only take neutral configu- 
rations. 

Introducing the charge density p(Q) = J2 a ,i q.a$(&—&a,i), the hamiltonian (2.12) 
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can be put in the form 



&a 2 v c{^a,ii &a,i 



(2.16) 



The term v c {Cl a ^, £l a ,i) which is subtracted corresponds to the term when, for 
the same particle, Q = Q' in the integral which should not be included in the 
hamiltonian. Using the explicit expression of the self energy (2.11), we notice 
that e a - qlv c ((l at i, £l a>i )/2 = -qlv (r a>i , r a>i )/2 is the self-energy of a particle 
in the flat geometry. 



Let us decompose the hamiltonian in spherical harmonics Yi m - Let 

oo I 

P(^) = E E pe m Y em (£l) . 
e=om=-e 



(2.17) 



Using (2.4) we have 



\Plr, 



oo I 

1=1 m=-l C \ l ^ L > a 



~Z!^ (r aii ,r aii )/2. 



(2.18) 



There are two problems to proceed to apply the sine-Gordon transformation 
to the grand partition function. First, the potential v c from equation (2.13), 
has no inverse. This can be seen in the expansion in spherical harmonics given 
by equation (2.4), with the term for £ = 0, V = 0. Second, to apply the sine- 
Gordon transformation one needs to consider all configurations of the system 
including the non-neutral ones. However, both problems are related and can 
be solved jointly. 



Let us write 5 



liuu 



,exp 



j3e{J2 a q a N a ) 2 j '4 in the partition 



function (2.15). Then we have S = lim^oo S e , where E e is a partition function 
not restricted to neutral configurations and with hamiltonian 



H e = H + e 



^q a N a 



/4- 



(2.19) 



This last term can be seen as a I = component to the interaction between 
pairs of particles. Indeed, since p o = R~ 2 E a ,i Q a Yoo(^a,i) = (47r)~ 1/2 i?~ 2 J2 a QaN a , 
we have 



H f = 7rR 



oo i I 12 
2 | \ ^ \ ^ \l J t,m\ 



(2.20) 



This solves also the second problem we faced before to apply the sine- Gordon 
transformation: the quadratic form in the hamiltonian has now a non- vanishing 
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£ = term and it is now invertible. At the end of the calculations we should 
take e — > oo. Thus, the £ = term of (2.20) goes to infinity. If we remember 
that the Coulomb potential is the inverse of the Laplacian operator, this is a 
reminder that !/[£(£ + 1)] diverges when £ = 0. This £ = term that we nat- 
urally introduce to take into account only the neutral configurations, has also 
modified the pair potential to render it invertible. Its inverse is the Laplacian, 
plus a £ = component proportional to 1/e. 

We now proceed as usual and perform the sine-Gordon transformation. In 



Refs. [8,9] we showed that the role of the flat self-energy t 



erm v 



in (2.20) is to 



regularize the ultraviolet divergence of the other terms. Therefore we will con- 



centrate our efforts only in the part H' t = irR A 
Performing the sine-Gordon transformation [ 



P00 e + ^£=1 Yl lrn =-t 

1], we have 



g _ jn^ n^ = -/^ An ex P [-g({c a })] 

Zo " / UT=o UL=-e d<t> lm exp[-S(0)] 



with the action 



(2.21) 



s({( a }) = -^T, 

£,m 



£m 



exp 



oo e 

£=0 m=-e 



:: R 2 dtt. (2.22) 



We have defined A lm = (3£(£+1)/(2tt) if £ > 1 and A m = (3e/(2n). (A tm ) repre- 
sents an operator which is the inverse of the modified Coulomb potential. It is 
proportional to the Laplacian plus a £ = component. = J2e, m <j>emYem($l) 
is the traditional auxiliary field introduced in the sine-Gordon transformation. 
The notation ::•••:: is defined as 



.. e -if3q a <P(Q) .. = e -if3q a <P(n) e f3q 2 a v°(r,r) 



(2.23) 



with r = (R, fl) , and can be understood as a sort of normal ordering operator, 
for details see [12]. 



The I — term in the action S is special, because at the end we will be 
interested in the limit e — > oo. Let us do the change of variable in the functional 



integral X = y / / 3/(27re)0 O o- The action is now 



^({Ca})= E ( t ) tm + — 

C>l,m 



(2.24) 



:: e 



-i/3q a 



Ei>i m <t>imY tm m+^/e/{2(3)X 



:: R 2 dQ. 
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2.3 Steepest descent 



We now consider the low-coupling or Debye-Hiickel regime, when (3e 2 <C 1, 
with e the elementary charge (all charges q a = z a e are supposed to be multiples 
of e, 4 e Z). In the low-coupling approximation, the action (2.24) is expanded 
to the quadratic order around the solution of the stationary action equation 
SS/5(p = 0, then the functional integral in equation (2.21) is Gaussian and 
can be performed exactly. The stationary action equation 5S/5(fi = reads, 

for e>i, 



P did , i\a , V" -a r fv fn\ ~ ipqa £*>i m <t>e m Y em (n)+^/e/(2(S)x 
— e(t+l)<i>i m +2^iPq a C a J Y em (Q)e 

and, for £ = 0, 

X + £i/??«C« J e -M^>-^ M ^ + -^ X } R^dn = (2.26) 



R 2 dVt = 
(2.25) 



Since the system is homogeneous and isotropic, we look for a constant solution, 
(pem = for £ > 1, which naturally satisfies (2.25). Then equation (2.26) 
reduces to 

— 1= + iR 2 £ q a Cae-^yffW> x = (2.27) 

Since we are interested in the limit e — > oo, the first term can be neglected 
and we have 

E fcCae-^v^)* = o (2.28) 

Notice that, if the fugacities are chosen to satisfy the pseudo neutrality con- 
dition J2 a <la(a = 0, then X = is a solution of (2.28). But X = X n = 



n27ry 2/(/3e 2 e) is also a solution with n£Z. Contrary to the usual situation 
in the infinite geometry or geometries with boundaries [8,9] when only the 
= stationary solution contributes, here all stationary solutions contribute 
to the functional integral (2.21), since when e — > oo, all X n — > 0. 

Let us proceed more generally, supposing that the fugacities are arbitrary and 
do not necessarily satisfy the pseudoneutrality condition. Let ip G K be the 
solution of 

E9«C«e-^ = 0, (2.29) 

a 

which is unique since the l.h.s. of (2.29) is a monotonous function of ip . Let us 
define (* = ( a e~P qa ^°. These "renormalized" fugacities satisfy the pseudoneu- 
trality condition J2 a C<?a — 0. Clearly the general solutions to equation (2.28) 
are 

X n = 2nn^^ + X , neZ (2.30) 
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with X = -iy/(2p/e)rl> . Now, we proceed to apply the steepest descent 
method to evaluate (2.21). Notice that it is valid to apply this method when 
f3e 2 <C 1 and ( a R 2 ^ 1- This last condition is necessary for the part of the 
integral depending on X. Taking into account all the stationary points X n we 
have 



1 



.-£/ n d^dxe-^-^) 2 ^ 2 ^ 

n&L J t>l,m 



Zy 

x exp 



(2.31) 



~ E [^ + l) + M) 2 ]::k:: 2 ^(l + ^) 2 )(^-^ 

£>l,m 



with k = ^j2ir(3J2 a Q 2 Ca the inverse Debye length. Performing the Gaussian 
integrals we find 



{nRf 



1-1/2 



n ( 1+ ^ri))n 



," 2 A2 



t>l,m 



e^„fi/( e ) (2.32) 



where \° k = — K 2 , K G M 2 , are the eigenvalues of the Laplacian operator in the 
flat space M 2 (these appear from the self-energy term v in the hamiltonian, 
for details see [8,9,12]), and V = 4nR 2 is the total area of the sphere. In the 
previous equation we have defined 



/(e) = (1 + {nRft)- 1 ' 2 J2 e-«(^ n -^) a 



which can be expressed in terms of the Jacobi function i? 3 (u|t) = J2nei> e 
as 

4iri 



/(e) = (1 + (KRfey^e^s ' 



ee 



/3e 2 e 



(2.33) 



iirrn 2 +2nui 



(2.34) 



To compute the limit of S e when e — > oo, it is useful to use the Jacobi imagi- 
nary transformation [13] to express /(e) as 



(3e 2 e 



\ 4tt(1 + M) 2 e) 



# 3 



i(3e 2 e 
Ait 



(2.35) 



The last term, the $3 function, has limit 1 when e — > 00. Finally, taking e — > 00 
we obtain the original partition function of the Coulomb system in the sphere 



-1/2 



yy 



(2.36) 



It should be clear to the reader that the above calculations are very general and 
can be easily adapted to other types of finite geometries without boundaries. 
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In that general case, the spherical harmonics and — l[l + 1)/R 2 are replaced, 
respectively, by the eigenfunctions and the eigenvalues of the Laplacian in the 
considered geometry. For a finite surface without boundaries, any constant 
function is an eigenfunction of the Laplacian with zero eigenvalue. This con- 
stant function plays the role of loo- Thus in an arbitrary finite surface without 
boundaries the grand partition function of the Coulomb gas is 




M e n,A n ^0 \ n I k 



-1/2 

e y £ a C. (2.37) 



where \ n are the eigenvalues of the Laplacian in the manifold where the system 
lives and V is the volume (area) of the manifold. The result is similar to the 
one for a geometry with boundaries and Dirichlet boundary conditions, found 
in Ref. [8], except that it appears an additional term K 2 V/([3e 2 ) which is the 
contribution of the zero eigenvalue. 

Notice that the role of the mean field ipo is to renormalize the fugacities, and 
the partition function is expressed in term of the fugacities (* which satisfy 
the pseudoneutrality condition. We recover the fact that the fugacities are 
not independent controlling variables, this is due to the global neutrality of 
the system, which we imposed from the start. The same situation arises in 
geometries with boundaries, for details see the appendix B of Ref. [8]. 



3 Thermodynamics of a Coulomb system on a sphere in the Debye— 
Hiickel regime 



3. 1 Grand potential 

In this section we consider the sphere geometry and we compute explicitly the 
partition function (2.36). The results of this section are specific to the sphere. 

The grand potential Q = — ksT \nE is given by 

-1>C (3-1) 

a 

The integral term in equation (3.1) comes from the eigenvalues \° k of the 
Laplacian in the flat space. The infinite product in equation (3.1) is divergent 
for large I and it should be cutoff to a maximum value iV of I. This ultraviolet 
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divergence is compensated by the one from the integral when K max — > oo. As 
explained and illustrated in several examples in Refs. [8,9], the cutoffs K mSuX 
and iV are proportional and the exact relation between them can be found by 
requiring that in the limit R — ■> oo we recover the bulk grand potential in the 
flat space which is known. In the integral we also need to introduce an infrared 
cutoff fcmin- It is explained in Refs. [8,9] that this cutoff is k m \ n = 2e~ c / L, where 
C is the Euler constant. The length L is the same from equation (2.2) which 
appears in the flat Coulomb potential v°. In the Debye-Hiickel regime it is 
understood that L is large, L — > oo. 

We can obtain an explicit expression for a certain regular izat ion of the infinite 
product appearing in (3.1). Let us consider the infinite product 

oo / 2 \ ^ 

- n (i + Kf+rj) < 3 ' 2 > 

which is convergent. Let us introduce the Barnes G function [14,15] 

G(z + 1) = ( 27 r)*/2 e -*(*+i)/2-c*»/2 J] ( 1 + 7 ) e~ z+z2 /W (3.3) 

with C the Euler constant. Let us write in equation (3.2) 

z 2 



+ 1) 

with 



1 + £ 



1/1,2 = (-1 ± Vl - 4^ 2 ) /2 (3.5) 

that satisfy £(£+l) + z 2 = — 1/1) — 1/ 2 ) . This allow us to express the product 
P(z) in terms of the Barnes G function, generalizing a standard procedure used 
to express infinite products in terms of Gamma functions [13], 

P( z ) = n (l + WT —-^ ) e- z 'l" = G(l - Vl )G{l - v 2 ) e-( 1+c > 2 (3.6) 
Then we have 



+ 1) 



TT (y 4. ^ ^ 2 ' +1 _ [^(1 ~ ~ ^ 2 )] 2 e- 2 (^ 2 

Ml 1 + ^TT)J " ra-^ra-^) (3 ' 7) 

Using r(| +x)r(| —x) = it/ cos(irx) and G(l + z) = T (z) G(z) we finally have 

ft i\ 4- *° ^ 2 ' +1 = vr[G(-, 1 )G(-, 2 )] 2 e^(^ 2 

M\ z 2 cosh([7rv / 4i^T]/2) 1 J 
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Putting z = kR and replacing (3.8) into the expression (3.1) for the grand 
potential yields 



= l - m W + 1 ln feiM^I (3 .9) 

2 /?e 2 2 ( Kj R)2 cos h([7r v /(2K J R)2-l]/2) 

+M) 2 f-i-^ + E7-inf^)-vEC« 

In the limit where the ultraviolet cutoffs i^ max and N go to infinity and re- 
placing k min = 2e~ c /L, we have 



^ = l ln il + l ln '[g(-^)G(-^)]- (3 . 10) 
2 /3e 2 2 cosh([ 7 r v / (2^R)2^1]/2) 



We see that if the ultraviolet cutoffs iV and K max are proportional the expres- 
sion for the grand potential is well defined. Furthermore, in the limit R — > oo 
we should recover the bulk value [81 



» K 2 [ KL 1 



EC (3.H) 



a 



for the grand potential. We can compute this thermodynamic limit using the 
known expansion of the Barnes G function for large argument [14,15,16] 

ln G(l + z) ~ z 2 - fj + Z - ln(2vr) - ^ + C'(-l) + 0(1/ z) . (3.12) 

1) is the derivative of the Riemann zeta function evaluated at —1). Using 
this expansion we find that it is necessary that K max = N/R to recover the 
correct bulk grand potential in the limit R — > oo. Finally the grand potential 
for finite R can be expressed as 



2 /3e 2 cosh([7r^(2 Kj R) 2 - l]/2) V L / 

-VEC (3-13) 

Using (3.12), we can find the finite-size expansion of the grand potential 

(3Q = (3Q b + l HkR) + 2C'(-1) - ^ - \ In ^ + 0{1/{kR)) (3.14) 
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with the bulk grand potential f2& given by equation (3.11). This expansion 
for the grand potential shows that there is a finite-size logarithmic correction 
(x/6)lni?, where x — 2 is the Euler characteristic of the sphere. Since the 
system has no boundary there is no surface (perimeter) tension term. 



3.2 Densities, pressure and the equation of state 

The density n a of the species a can be obtained using the standard thermo- 
dynamic relation n a = ( a d\nE/d( a . When computing this derivative one has 
to take into account the fact that ipo is a function of ( a implicitly given by 
equation (2.29). After some long but straightforward calculation we finally 



where ip(x) = d In T(x) / dx is the psi function. This relation allows us to express 
the derivative of the Barnes G function in terms of the psi function ip(x). 
Also we used some known identities [10] satisfied by the psi function, such as 
ipiy) + ip{—u — 1) = 2ip(u + 1) — \v{y + 1)] _1 + 7r cot(7rz/). One can easily verify 
that the system is neutral J2 a Qan a = 0. The total density n = J2 a n a is 



find 




(3.15) 



We have used the relation [15] 




(3.16) 




(3.17) 



The finite-size expansion of the densities, when R — > oo, reads 




(3.18) 



where the bulk density is 




(3.19) 
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And for the total density 



with n b = J2 a n a- Thus the total number of particles N = N b + 1/6, where N b 
is the bulk number of particles. The 1/6 finite-size correction to the number of 
particles is a consequence of the (1/3) lni? correction to the grand potential. 

The pressure p = —dVt/dV is given by 

fo = EC- ^ + ^^(l + ^i) + fcot(^ 1 ) + In ^p) • (3.21) 

Using equation (3.15), and neglecting terms of higher order than j3e 2 , we find 
the equation of state 

0P = X>«(l-^) • (3-22) 

We remind the reader that, from a scale invariance analysis, one can show 
that this equation of state is actually valid in the whole range of stability of 
the system of point particles, both in the flat geometry and in the sphere. 



3.3 Internal energy 



The excess internal energy is U cxc = (d(/3Q)/d/3)^ a y. Using (3.13), we find 

(3U exc = -(kR) 2 (V(l + vi) + | cot(7rz/ 1 ) + C + In ^ _ I ( 3 . 23 ) 
We recall that v x = (-1 + y/l - (2kR) 2 ) /2. 

We can derive this result (3.23) using a more traditional approach to the 
Debye-Hiickel theory. Consider a unit pseudocharge (charge plus neutralizing 
background) located at the north pole. Following the usual formulation of 
Debye-Hiickel theory, this pseudocharge is screened by a polarization cloud 
created by the plasma. This cloud has a charge density given by 



~-K 2 D K{0,<p;Q,G) + K 2 D (K) (3.24) 

where K(0, (p; 0, 0) is the (mean field) electric potential created at (8, ip) by the 
pseudocharge and its polarization cloud. Also (K) = J K{6' ', <//; 0, 0) dQ! / (47r) 
is the average of K over the sphere. To understand the second term in (3.24) 
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recall that we are dealing with pseudocharges (charges plus neutralizing back- 
ground), and the polarization cloud is also made of pseudocharges. In the 
second line we linearized the exponentials as it is usually done in the Debye- 
Hiickel theory and defined the inverse Debye length kd — \J^f3 J2 a n a q^. At 
the Debye-Hiickel level of approximation k d ~ k. The potential K satisfies 
the modified Poisson equation 

AK = -27T (5 - — + Ppol (9, ^)) (3.25) 

and using (3.24) we arrive at the modified Debye-Huckel equation for the 
sphere geometry 

AK-k 2 K + k 2 (K) = - 2% ( 6 -^]^) ( 3 - 26 ) 

As in the case of the modified Poisson equation (2.3), we can notice that the 
solution of the modified Debye-Hiickel equation (3.26) is determined up to an 
arbitrary additive constant. Indeed this can be clearly seen if we look for a 
solution of equation (3.26) as an expansion in spherical harmonics 



K(6, V ; 0,0) = Y, E ,,, lU , B J 40,0my) (3-27) 

e=lm=-£ l \ l + l > + \ Ktt ) 

+& o y O o(o,o)y O o(0,<p)- 

The constant term & , for I — 0, cannot be determined from the differential 
equation (3.26). Thus the average (K) = 6 /(47r) can be chosen arbitrarily. 
However as we will see below this constant term will not be needed in the 
following. 

The excess internal energy ?7 exc can be computed from the Debye-Hiickel po- 
tential K as the potential energy of the pseudocharge in the north pole, in the 
potential field K 



U exa =——Y t ^ a n a (3-28) 



hm 



J K(6, ip; 0, 0)(5(9, if; 0, 0) - 1 / (Ait R 2 )) R 2 dtt + In % - 



Notice that we subtract the self-energy — ln(r/L) corresponding to a flat ge- 
ometry, in accordance to the prescription (2.10), to obtain a finite result. This 
kind of prescription has also been used in Ref. [17] for the derivation of the 
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pressure of the Coulomb gas using the Maxwell stress tensor. We have 
4ttR 2 



U e xc 



J2qln a lim (k(6,<p; 0, 0) - (K) + In £j (3.29) 



Thus we only need K — (K), the constant term (K) is irrelevant for the 
calculation of the internal energy. An explicit expression for K can be obtained 
by noticing that the Yukawa potential Y, which is the solution of 



AY - k 2 Y = -2rx5 , 



(3.30) 



is also a particular solution of equation (3.26) with (Y) = 1/{2k 2 R 2 ). There- 
fore, Y has an expansion in spherical harmonics of the form (3.27) with 
b = 2-k/(k, 2 R 2 ). So, the difference between K and Y is a constant equal to 
-1/[2(kR) 2 ] + (K). On the other hand equation (3.30) can be solved directly 
since it reduces to a Legendre equation in the variable cos 9. The solution 
is [18] 



Y{9) 



71 



2 sm(i/i7r) 



cos 9) 



(3.31) 



with P Ul a Legendre function. Finally we find the Debye-Hiickel potential for 
a pseudocharge located at the north pole 



-P Ul (-cos9)- ' 



2sin(z/ 1 7r) " 1V ' 2(kR) 2 

As 9 — > 0, the Legendre function has the behavior [18] 



(3.32) 



P V1 (- cos( 



2 sin(i/i7r) 



7T 



9 7T 

In sin — h C + + v x ) H — cot(i/i7r) 
2 2 



0(1). 

(3.33) 

Using this asymptotic behavior into equation (3.29) allow us to retrieve the 
result (3.23) for the internal energy. 



4 Finite-size corrections for Coulomb systems on finite surfaces 
without boundaries 



In this section we consider a Coulomb gas confined in an arbitrary finite sur- 
face without boundaries. We will compute the finite-size expansion of the 
grand potential. Let R be the square root of the total area V of the surface. 
Following [9], we introduce the zeta function of the Laplacian in the geometry 
considered, 

oo 

Z(s,a) = Y J {a-\ k )- s (4.1) 
k=i 
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where are the eigenvalues of the Laplacian. Notice that we omit the van- 
ishing eigenvalue A = in the definition (4.1). 



In Ref. [9], this zeta function is related to grand potential. We can directly 
transpose the calculations of Ref. [9] to the present case, taking special care 
of the additional contribution to the grand potential due to the vanishing 
eigenvalue. Then, we obtain for the grand potential the relation 

iM = \ [Z'(0, 0) - Z'(0, k 2 )] + ^ In ^ + 1 In - V £ Q- (4.2) 

The prime means differentiation with respect to first variable of the zeta func- 
tion. The finite-size expansion of the zeta functions involved in (4.2) can 
be obtained from the known small-argument expansion of the heat kernel, 
9(t) = E^ =0 e Ufc , of the Laplacian [19,20] which reads, 

0(t) = ^ + | + o(tV 2 ) (4.3) 

with x the Euler characteristic of the surface. This is explained in detail in 
Ref. [9]. Here we only need to take special care of the fact that we omitted 
the zero eigenvalue in the zeta function, which is equivalent to subtract 1 to 
the heat kernel. Thus, following [9], we find 



1 r Z'(0,0)-Z'(0,K 2 )] =|l(i-in«)y+ x 



- - 1 

6 



MkK) + 0(1) (4.4) 



Replacing into (4.2) we finally find the finite-size expansion of the grand po- 
tential 

(3n = (3n b + ^\n(KR) +0(1) (4.5) 

with the bulk grand potential given by (3.11). Notice the existence of a 
logarithmic finite-size correction (x/Q)lnR, which is universal, i. e. indepen- 
dent of details of the microscopic constitution of the system. Actually this 
universal finite-size correction seems to exist even beyond the low-coupling 
regime considered here. In particular it has been shown to exist for the sphere 
geometry for the one-component plasma [21] and the two-component plasma 
both in its charge symmetric [22] and charge asymmetric version [23], in the 
whole range of stability of the system of point particles. 



5 Summary and perspectives 



We have showed how to build the Debye-Huckel theory for two-dimensional 
Coulomb systems confined in a finite surface without boundary. In particu- 
lar we showed how to perform the sine-Gordon transformation for this kind 
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geometry and how to map the statistical mechanics problem into a field the- 
ory on the finite surface. This could have further applications, for instance, 
to compute higher order corrections in f3e 2 to the grand potential and other 
thermodynamic quantities. 

For the case of the sphere geometry we explicitly computed the grand poten- 
tial and other thermodynamic quantities. For a general geometry, we showed 
the existence of a universal logarithmic finite-size correction for the grand 
potential. 
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